Soft Open Charm Production in Heavy-Ion Collisions. 
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Effects of strong longitudinaf cofor electric fields (SCF) on the open charm production in nucleus- 
nucleus {A + A) collisions at 200A GeV are investigated within the framework of the HIJING/BB 
v2.0 model. A three fold increase of the effective string tension due to in medium effects in ^4 + j4 
collisions, results in a sizeable (~ 60-70 %) enhancement of the total charm production cross sections 
{(J^^). The nuclear modification factors show a suppression at moderate transverse momentum (pt) 
consistent with RHIC data. At Large Hadron Collider energies the model predicts an increase of 
by approximately an order of magnitude. 
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The phase transition from hadronic degrees of free- 
dom to partonic degrees of freedom in ultra-relativistic 
nuclear collisions is a central focus of experiments at the 
Relativistic Heavy Ion Collider (RHIC). Heavy- flavor 
quarks are an ideal probe to study early dynamics in 
these nuclear collisions. Several theoretical studies pre- 
dict [J [2, 9 ^ substantial enhancement of open charm 
production associated to plasma formation of the de- 
confined parton matter relative to the case of a purely 
hadronic scenario without plasma formation. A recent 
analysis shows that the dynamics of heavy-quarks at 
RHIC are dominated by partonic or "pre-hadronic" in- 
teractions in the strongly coupled plasma (sQGP) stage 
and can neither be modeled by "hadronic interactions" 
nor described appropriately by color screening alone [1] . 
Therefore, these quarks are key observables in the study 
of thermalization of the initially created hot nuclear 
matter Q. 

A review of heavy-flavor production in heavy-ion col- 
lisions has been recently published Direct recon- 
structed mesons via hadronic channel (D° Kit) 
in d + kvL 01, Cu + Cu and Au Au @ colli- 
sions have been measured. Due to the difficulty to re- 
construct D-mesons hadronic decay vertex, both STAR 
and PHENIX have studied open charm indirectly via 
semileptonic decay to non-photonic electrons (NPE) or 
muons [1 i, E, 111 [II, [S B [II- Theory predict 
that charm quarks are produced by initial gluon fu- 
sion [l^ and their production rates are expected to 
be well described by perturbative Quantum Cromody- 
namics (pQCD) at Fixed Order plus Next to-leading 
Logarithms (FONLL) [13 . Total charm cross sections 
reported by both experiments are however only com- 
patible with the upper limit of the FONLL predictions. 
In addition, the data indicate a suppression as large as 
that of hght quarks [I3|,[IBl while, due to their large 
mass and to the dead cone effect charm quarks are pre- 
dicted to loose less energy than light quarks by gluon 
radiation in the medium |18l |. 

Recent model calculations based on in-medium charm 
resonances/diffusion or collisional dissociation [I^, [2^, 



l2l| . radiative energy loss via few hard scatterings 123 
or radiative energy loss via multiple soft collisions [23|, 
have been applied to describe the non-photonic elec- 
trons (NPE) spectra. They all predict less suppression 
than that observed in experiments. On the other hand, 
a good description of the nuclear modification factor 
(NMF), RJ2^(pt), was obtained in non-perturbative 
time-dependent heavy-quark diffusion in the Quark- 
Gluon Plasma (QGP) [3- 

In previous papers [2^,[2^ we have shown that the dy- 
namics of strangeness production deviates considerably 
from calculations based on Schwinger-like estimates for 
homogeneous and constant color fields [131 and point 
to the contribution of fluctuations of transient strong 
color fields (SCF). These fields are similar with those 
which could appear in a ^^glasma" [2^ at initial stage 
of the collisions. In a scenario with QGP phase transi- 
tions the typical field strength of SCF at RHIC energies 
was predicted to be about 5-12 GeV/fm 12911. Recently 
Schwinger mechanism has been revisited j30'| and pair 
production in time-dependent electric fields has been 
studied [31j]. It is concluded that particles with large 
momentum are likely to have been created earlier than 
particle with small momentum and for very short tem- 
poral widths (Ar « lOtc, where the Compton time 
tc — 1/mc) the Schwinger formula strongly underes- 
timates the reachable particle number density. 

In this paper we extend our study in the framework of 
HIJING/BB v2.0 model [2^ to open charm productions. 
We explore dynamical effects associated with long range 
coherent fields (i.e strong color fields, SCF), including 
baryon junctions and loops |25| . with emphasis on the 
novel open charm observables measured at RHIC in p+p 
and heavy-ion collisions. Using this model we analyze 
the enhancement of total charm production at 200A 
GeV energy. 

For a uniform chromoelectric flux tube with field {E) 
the pair production rate [sol [s^ . [33| per unit volume for 
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a heavy quark is given by: 



k2 / 7rm^ 



(1) 



where for Q = c or 6 , mp = 1.27, or 4.16 GeV 
(with ±1% uncertainty [Si) Note that k ^ \eE\eff = 
y/ C2 {A) IC2 (F) kq is the effective string tension in 
terms of the vacuum string tension kq ~ 1 GeV/fm and 
C2{A), C2(F) are the second order Casimir operators 
(see Ref. [33] )■ In a nuclear collisions, the local lon- 
gitudinal field strength increases with the square root 
of the number of color exchanges proportional to the 
number of binary collisions per unit area {n{x±,b) cx 
^jTAA{x±^h)), where for a given impact parameter 
b and transverse coordinate xj^, TaA cx A^^^ is the 
Glauber A+A binary collision distribution. Therefore, 
the effective string tension k c>c A^^^ . 

A measurable rate for spontaneous pair produc- 
tion requires '''' strong chromo electric fields^^ , such that 
k/ttiq > 1 at least some of the time. On the aver- 
age, longitudinal electric field "string" models predict 
for heavier flavor a very suppressed production rate per 
unit volume 7q via the well known Schwinger formula 
[131, since 



exp 



7r(m2 - ml) 



« 1 (2) 



for Q = c and q — u,d. For a color rope on the other 
hand, if the average string tension value (< k >) in- 
creases from 1.0 GeV/fm to 3.0 GeV/fm, the rate F for 
charm pairs to tunnel through the longitudinal field in- 
creases from « 1.4 • 10~^^ to « 3.5 • 10~^ fm""*, and this 
can lead to a net "soft" tunneling production compara- 
ble to the initial "hard" FONLL pQCD production. 

The conventional hard pQCD mechanism, mainly 
gluon fusion Jj , is calculated via the PYTHIA subrou- 
tines in HIJING/BB v2.0. The advantage of HIJING over 
PYTHIA is the ability to include novel SCF color rope 
effects that arise from longitudinal fields amplified by 
the random walk in color space of the high x valence par- 
tons in j4-|-^ collisions. This random walk could induce 
a very broad fluctuation spectrum of the effective string 
tension. Thus, if the average or mean < k >= nno, 
then the typical fluctuation is of order 1/Vn which is 
large because n « 6 for Au nuclei. A Poisson fluctua- 
tion of effective k about the mean, gives a strong bias 
toward less probable but larger than the average value 
< K >. This is amplified for heavy quarks. Here we 
do not investigate in details such fluctuations, but we 
will estimate the effects of a larger effective value k > 
3 GeV/fm on the enhancement of . 

Both STAR and PHENIX experiments have mea- 
sured charm production cross sections in several colli- 
sion systems. Figure [T] shows the measured total charm 
production cross sections at mid-rapidity, da^^/dy 



(left panel) and in all phase space, (right panel). 
The data from both experiments seems to indicate a 
scaling with number of binary collisions (Nbin), as ex- 
pected because of the high mass of charm pairs pro- 
duced in initial nuclcon-nuclcon collisions [16|. How- 
ever, there is still an unresolved discrepancy of the order 
of a factor of two between STAR and PHENIX data. 
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FIG. 1: (Color online) Comparison of HIJING/BB v2.0 pre- 
dictions for mid-rapidity (left panel) and all phase space 
(right panel) charm cross sections per nucleon-nucleon col- 
lisions as a function of A'^bin in {(I)A-\-A collisions. The sym- 
bols are the results with (filled squares) and without (open 
crosses) SCF efi'ects. Both include quenching and shadow- 
ing (ys) efi'ects. The open triangles are the results with SCF 
effects but no shadowing (ns) . The values of FONLL predic- 
tions are shown as a dotted line. The band at the left mark 
the FONLL uncertainties [13] ■ The data are from STAR 
(stars) [3|,[l,[il,[nil and PHENIX (solid circles) [il,[Il. 
Statistical and systematical error bars are shown. 



The predictions of HIJING/BB v2.0 model without 
SCF (open crosses) and including SCF effects (filled 
squares) are shown in the figure. For completeness 
the results with SCF but no gluon shadowing effects 
(open triangles) are also included. However, in this sce- 
nario multiplicities at mid-rapidity are strongly over- 
estimated [35| . The main parameters used in the cal- 
culations are given in Table II of reference [2^, and 
corresponds to strengths of strong color (electric) field 
dependent on collision system (k = 1.5; 2.0; 3.0 GeV/fm 
for p + d~\- Au, and A + A collisions respectively) . In 
our calculations we estimate the total open charm pro- 
duction (c -I- c) cross section considering the 12 lightest 
I)-mesons (0°, 5°, 0°*, 5°*, D+, D+, p+*, D+*, D^, 
Dg, D*, D*), and the hyperons Ac and Ac. The contri- 
bution of higher mass charm hyperons is negligible. For 
calculations which take into consideration SCF effects 
(filled squares) we obtain an increase of 60 — 70 % in 
comparison with a scenario without SCF effects (open 
crosses). These results describe well the PHENIX data 
within statistical and systematical errors and are close 
to the upper limit of uncertainty band of the pQCD 
FNOLL predictions [13] ■ Our calculations also show 
that the scaling with Nbin is only approximately sat- 
isfied, the reason being an interplay between the mass 
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dependent SCF and shadowing effects, which act in op- 
posite directions. In fact, we calculate that only 60 % of 
total open charm production (c + c) conies from par- 
tons embedded within the target and projectile. 

The study of open charm production in rf + Au col- 
lisions allow to separate "cold nuclear matter" effects. 
The initial production of c c pairs by gluon fusion might 
be suppressed due to gluon shadowing. We recall that 
shadowing is a depletion of the low-momentum parton 
distribution in a nuclcon embedded in a nucleus com- 
pared to a free nucleon; this leads to a lowering in the 
(scaled) c + c production relative to p + p collisions. 
The shadowing in the regular HIJING parameteriza- 
tion [s^l implemented also in our model seems to be 
too strong [37] . There is a considerable uncertainty (up 
to a factor of 3) in the amount of shadowing predicted 
at RHIC energies by the different models with HIJING 
predicting the strongest effect [s^. This could explain 
why the results for scaled cross sections in d -I- Au colli- 
sions are smaller than those obtained for p + p collisions 
(see Fig. [T] left panel) . 

We study if we can find scenarios that would give 
larger enhancement of total cross sections for open 
charm production, than those reported in Fig. [T] (filled 
squares), and that would be consistent with the STAR 
data. The random walk in color space of heavy quark 
could induce a broad spectrum of the effective string 
tension. Therefore, we study the effects of a further 
increase of mean value of the string tension from 3.0 
GeV/fm to 5.0 GcV/fm on c c pair production. This 
results in only a modest increase of scaled cross sec- 
tions, af^ = (T^^/Nbin, by approximately 20% for 
central collisions. For values between 5-10 GeV/fm a 
saturation sets in, as an effect of energy and momen- 
tum conservation constraints. In our model the mul- 
tiplicative factor that accounts for next to leading or- 
der corrections [s^ in calculations of hard or semihard 
parton scattering processes via pQCD is set to K = 2. 
Increasing this factor to K = 3.5, as suggested in refer- 
ence [3^ results in an increase of cr^Jf^ by approximately 
70% in central Au + Au collisions, but also overpredicts 
by 40 % the total charged particle production at mid- 
rapidity (Nch (y=0)). Therefore, we conclude that the 
large charm cross sections obtained by the STAR collab- 
oration cannot be explained within our phenomenology. 

The £) '^-mesons spectra are sensitive to the dynamics 
of produced charm particles. In Fig. [2] we present the 
calculated D°-mesons spectra for systems where data 
are available jij,iS],[S]. In all cases, the calculated yield 
is much smaller than the STAR data, consistent with 
the results shown in Fig.[TJ The calculated spectra show 
little shoulder at low pt indicating small radial flow of 
D"-mesons consistent with the results of STAR [llj . 

Figure [3] shows our predictions for the Nuclear Modi- 
fication factor (NMF), Raa(pt) for and 7r° mesons. 
Data (filled symbols) are NMF for non photonic elec- 
trons, R^']^(pt) [i3|j[IB|- The data for tt" meson (open 
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FIG. 2: (Color online) Comparison of HIJING/BB v2.0 pre- 
dictions for Pt distribution of invariant yields for recon- 
structed in minimum-bias {d)A+A collisions. For clarity 
the results for Cu -I- Cu (preliminary) and Au -I- Au are mul- 
tiplied by 10 and 10'^ respectively. The data are from STAR 
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symbols) are from reference [40|. Note, that non pho- 
tonic electrons include also electrons from bottom (b) 
production (B IX) and the yields of mesons could 
be affected by the decay (B D). For central (0-10 %) 
Au + Au collisions we calculate a scaled total cross sec- 
tion for bottom production with (without) SCF of (J^f^ 
= 17.8 ^b (cr^^^ = 0.86 /Ltb). These values are few orders 
of magnitude lower than af^ and this contribution is 
estimated to be negligible for px < 6.0 GeV/c. 
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FIG. 3: (Color online) Comparison of HIJING/BB v2.0 pre- 
dictions of nuclear modification factor Raa(pt) for D*^ and 
tt'' mesons in central (0-12 %) Au -I- Au collisions. Data 
from STAR (stars) ^ and PHENIX (circles) [3 are NMF 
for non photonic electrons, R^'^(pt). The data for tv" me- 
son are from PHENIX [Z^]. Error bars show the statistical 
and systematic uncertainties. 



In our calculations for low px (0 < px < 2.5GeV/c), 
non-perturbative production mechanism via SCF re- 
sults in a split between 1)° and 7r° mesons. The charged 
and TT^ mesons are suppressed due to conservation of 
energy [2^ . The yields of the I?'^-mesons are enhanced 
due to an increase of cc pair production rate (see Eq. 
1). In central (0 — 10%) Au -f- Au collisions a suppres- 
sion at moderate px (4 < px < 6 GeV/c) as large 



4 



as that of light quarks is observed in contrast to previ- 
ous theoretical studies [H , [H [H, [H, EH . Our model 
predicts a suppression consistent with the data. We 
can interpret this results as experimental evidence for 
"m-merfmm mass modification^^ of charm quark, due to 
possible induced chiral symmetry restoration |42l |. An 
in-medium mass modification has also been predicted 
near the phase boundary (i.e. at lower energy) in [43l |. 
In contrast statistical hadronization model |44l | predicts 
no medium effects at top RHIC energy. 

We performed calculations at the much higher Large 
Hadrons Collider (LHC) energy using parameters from 
reference i.e n = 2.0; 5.0 GeV/fm for p + p and cen- 
tral (0-10 %) Pb -|- Pb coUisions respectively. The pre- 
dicted charm production cross section is approximately 
an order of magnitude larger than at RHIC energy. We 
obtain = 6.4 mb in p-|-p collisions and a (scaled) 
cross section afl^ — 2.8 mb for central Pb + Pb colli- 
sions (Nbin = 960 and Nch (y=0) = 2500). This indi- 
cates a clear violation of scaling with Nbin at the LHC. 
These values increase by a factor of 2 to 3 if the effects 
of shadowing are not included ( Nch (y=0) ~ 5000 and 
aff w 8.4 mb). 

In summary, we studied the influence of possible 
strong homogenous constant color fields in open charm 
production in heavy-ion collisions by varying the ef- 
fective string tension that control QQ pair creation 



rates. This is equivalent with assuming an in medium 
mass modification of charm quark. We show that this 
approach is an important dynamical mechanism that 
can explain the observed _D-mesons enhancement pro- 
duction observed by the PHENIX experiments. Our 
model is based on the time-independent color field 
while in reality the production of QQ pairs is a far- 
from-equilibrium, time-dependent phenomenon. Thus 
to achieve more quantitative conclusions, such mecha- 
nisms [3l| should be considered in future generation of 
Monte Carlo codes. 

The large cross sections reported by the STAR collab- 
oration remain unexplained within our study. Solving 
the discrepancy between the measurements is impor- 
tant, since confirmation of the STAR results may in- 
dicate the importance of other dynamical mechanisms 
such as pre-equilibrium production from secondary par- 
ton cascades lll|, or hot-glue scenario Q. 
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